Statistics of current fluctuations at non-zero frequencies 
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We formulate a general approach which describes statistics of current fluctuations in mesoscopic 
coherent conductors at arbitrary frequencies and in the presence of interactions. Applying this 
approach to the non-interacting case, we analyze frequency dispersion of the third cumulant of the 
current operator S3 at frequencies well below both the inverse charge relaxation time and the inverse 
electron dwell time. This dispersion turns out to be important in the frequency range comparable 
to applied voltages. For comparatively transparent conductors it may lead to the sign change of S3. 



Investigations of current fluctuations in mesoscopic 
conductors can provide a great deal of information about 
properties of such systems. During last years much at- 
tention has been devoted to shot noisei. Experimental 
and theoretical studies of the second moment of the cur- 
rent operator describing shot noise revealed a rich vari- 
ety of properties caused by an interplay between scat- 
tering, quantum coherence and charge discreteness. Fur- 
thermore, such studies can substantially deepen our un- 
derstanding of the role of electron-electron interactions in 
mesoscopic transport because shot noise and interaction 
effects are known to be closely related 2 * 3 -. 

One can also go beyond the second moment and study 
higher order correlators of the current operator thereby 
extending the amount of information already obtained 
from investigations of electron transport and shot noise. 
Recently the first experimental study of the third current 
cumulant in mesoscopic tunnel junctions was reported- . 

A theoretical framework which enables one to analyze 
statistics of charge transfer in mesoscopic conductors was 
developed in Ref. 0. This theory of full counting statis- 
tics (FCS) allows to evaluate any cumulant of the cur- 
rent operator in the absence of interactions and in the 
zero frequency limit. Under these conditions higher or- 
der current cumulants were investigated by a number of 
authorsSiSiS. In order to include interactions and to 
analyze frequency dispersion of current fluctuations it is 
necessary to go beyond the FCS theory and to develop a 
more general real time path integral techniqu e 2 ' 10 . 

The goal of the present paper is to address statistics 
of current fluctuations at non-zero frequencies. We first 
present a general and formally exact expression for the 
real time effective action of a coherent conductor de- 
scribed by an arbitrary energy independent scattering 
matrix. This expression enables one to fully describe in- 
teraction effects in such type of conductors. We will then 
demonstrate that in a non-interacting case this effective 
action provides a direct generalization of the FCS gener- 
ating function* to non-zero frequencies. With the aid of 
our technique we will analyze the frequency dispersion of 
the third cumulant of the current operator in mesoscopic 
coherent conductors. 

It is worthwhile to point out that the frequency de- 
pendence of current correlators can be caused by various 



reasons. One of them is the effect of an external elec- 
tromagnetic environment which is important for quanti- 
tative interpretation of the experimentally detected be- 
havior of higher cumulantaii. Another source of the fre- 
quency dispersion is the internal dynamics of a quantum 
scatterer. Here the important time scales are the corre- 
sponding RC -time trc and the electron dwell time td 
inside the conductor. The latter scale was recently taken 
into consideration in the analysis of the second^ and the 
third 13 current cumulants for chaotic quantum dots. 

In this paper we will address current fluctuations at 
frequencies not directly related to any of such scales. We 
will demonstrate that apart from the above mechanisms 
there exists an additional - intrinsic - frequency disper- 
sion of the current correlators at the scale set by the 
voltage drop V across the conductor. Since V can vary 
in a wide range, this dependence is in general important 
and should be taken into account while interpreting the 
experimental results. In particular, in the absence of in- 
teractions the third cumulant of the current operator S3 
is fully determined by the two parameters: 
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where T n represents the transmission of the n-th con- 
ducting channel of our system. The cumulant S3 can be 
expressed in the following general form 



S 3 = 09 - 2 1 F)e 2 I. 



(2) 



Here I is the average current through the conductor, — e 
stands for the electron charge, and F is a universal func- 
tion of frequencies, voltage V and temperature T to be 
evaluated below. According to Eq. (J2J the frequency dis- 
persion of S 3 originates only from the term proportional 
to the parameter 7 while the /3-term is dispersionless. 

Though negligible for tunnel junctions 7 — > 0, the fre- 
quency dispersion of S3 may become important in other 
situations. For instance, at T — > one finds F — ► 1 at fre- 
quencies much smaller than eV while in the opposite high 
frequency limit one gets F — and, hence, S3 = (3e 2 I 
in the latter limit. To give some numbers, for an impor- 
tant case of diffusive conductors one has (3 = 1/3 and 
/? — 27 = 1/15, i.e. in this case the quantity S3 changes 



2 



by the factor 5 depending on whether relevant frequen- 
cies are below or above eV. For conductors with j3 < 27 
even the sign of S3 will differ in these two limits. 

In our analysis we will use the real time path inte- 
gral formalism developed for the systems of interacting 
fermions 14 . After the standard Hubbard-Stratonovich 
decoupling of the interaction term in the Hamiltonian 
one can exactly integrate out fermions and arrive at the 
effective action S which depends on the fluctuating fields 
^1,2 Let us define 



= Tr We~ 1 ^ dt ' H ^ t ') p^fe^o dt'H 2 (t') 



(3) 



with the trace taken over the fermionic variables. Here 
Po is the initial N— particle density matrix of electrons, 

H 1} 2 = £ J d 3 r¥ a (r)H h2 (t)^ a (r), 



H 1 , 2 {t) = -— + U{r)-eV ia {t,r) (4) 
2m 

are the effective Hamiltonians on the forward and back- 
ward parts of the Keldysh contour, U(r) describes the 
static potential, T and T are respectively the forward 
and backward time ordering operators. Integrating out 
fermions in Eq. ©, we obtain iS = iSg + iS cm where 



iS a = 2Trln[l + (u 2 u\ — l)p ], 



(5) 



wi,2(0 — ^~exp ^— i f* dt' Hx^it'yj are the evolution op- 
erators pertaining to the Hamiltonians @, po stands for 



the initial single-particle density matrix and the term S cm 
accounts for the electromagnetic contribution which may 
also include the effect of an external circuit. 

In order to evaluate the evolution operators 1*1,2 it is 
necessary to specify the model of a mesoscopic conduc- 
tor. Here we will adopt the standard model of a (com- 
paratively short) coherent conductor placed in-between 
two bulk metallic reservoirs. The electron dwell time tjj 
is supposed to be shorter than any relevant time scale in 
our problem. Energy and phase relaxation times are, on 
the contrary, assumed to be long, i.e. inelastic relaxation 
is allowed in the reservoirs but not inside the conductor. 
Under these assumptions electron transport through the 
conductor can be described by the energy independent 
scattering matrix 



S 



r t' 
t f> 



(6) 



and the effective action (jSJ) can be expressed via the 
fluctuating phase fields <p\ l2 which are in turn related 
to the jumps of the fields Vi }2 across the scatterer as 
^1,2 = g(Vli,2 — Vri,2)j where Vl,r are fluctuating in 
time but constant in space fields in the left and right 
reservoirs. We note that in this case the right hand side 
of Eq. © differs from the FCS generating functional 
introduced in Ref. only by a gauge transformation. 

Within the above mo del th e evolution operators u\p 
were evaluated in Refs. I2ll0i Combining these expres- 
sions with Eq. JSJ) after some algebra we find 



iS = 2Trln <^ 16(x - y) + 6(t - x)6(x) 



ftf( e ^ (*) - l) 
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Here we introduced ip + = (tpi + tp 2 )/2, ip = ipi — ip 2 and 



Po(x) 



dE 



■iEx 



2tt 1 + e E ' T ' 



(8) 



Taking the trace in Eq. (JJJ implies convolution with 
respect to internal time variables (x and y). The total 
time span is denoted by t. 

Eq. defines a formally exact effective action for 
a coherent conductor described by an arbitrary energy 
independent scattering matrix @. This expression al- 
lows to fully determine statistics of current fluctuations 
at arbitrary frequencies and in the presence of interac- 
tions. In the case of equilibrium fluctuations the formula 
Q represents a real time analogue of the effective action 
derived in Refs. ITiillTTl within the Matsubara technique. 
A formula similar to Eq. was also presented recently 



in Ref. jig. In addition we note that, provided the field 
ip~ does not depend on time, Eq. Q coincides with the 
generating function considered, e.g., in the problem of 
adiabatic pumping through mesoscopic conductors 1 ^. 

Let us illustrate the relation between the effective ac- 
tion and the FCS generating function^. For this pur- 
pose we (i) disregard interactions and set tp + {t) = eVt 
and (m) suppress fluctuations of ip~ and set ip~ — const. 
After these simplifications from Eq. J7J we obtain 



"iSfcs = — Tr / dE In 



l + iHn L (E)(l~n R (E)) 



x (e 1 ^ - 1) + #tn R (E)(l - n L (E))(e-^~ - 1)J , (9) 

where n LM (E) = 1/[1 + exp[(E ± eV/2)/T]]. Eq. © is 
just the FCS generating function- which can be used to 
recover all cumulants of the current operator in the zero 



3 



frequency limit and in the absence of interactions. 

The expression (jTJ) is more general since it includes all 
possible fluctuations of the phase fields ip^ . For instance, 
if one allows for temporal variations of ip~ , with the aid of 
Eq. (JJJ one can easily describe the frequency dispersion 
of the current correlators. Below we will illustrate this 
point by directly evaluating the third cumulant of the 
current operator. 

Within our model the current operator can be defined 
in a standard manner 



I(t) = - Jt (N L (t)-N R (t)), 



(10) 



where Nl(r) is the total number of electrons in the left 
(right) lead. Combining this definition with Eq. © for 
the expectation value one finds 

(I(t)) = -e [ V^SSl^ySfp-We* 3 ^. 



Similarly one can define higher moments of the current 
operator. Here we consider the following correlation func- 
tions: S 2 = (J(ti)/(t a ) + /(t a )/(ti))/2 and 

S 3 = ^{(7(t 1 )(TJ(t 2 )/(t 3 ))> + ((t/(t 2 )/(t 3 ))/(t 1 )) 

+ (i(t 2 )(Ti( tl )i(t 3 ))) + {(fi( tl )i(t 3 ))i(t 2 )) 
+ (i(t 3 )(r/(i 1 )7(t 2 ))) + ((r/(ti)/(t 2 ))/(t 3 )> 
+ (Ti( tl )i(t 2 )i(t 3 )) + (Ti(ti)i(t2)i(t 3 )}} . (ii) 

The correlation function S 2 is important because the 
symmetric combination of voltages V + can be viewed 
as a classical, measurable, voltage 15 . The noise is de- 
duced from the measurable product V + (ti)V + (t2), which 
is related to the symmetric correlator 52- Similarly the 
measured product V + (ti)V + (t2)V + (t 3 ) is related to the 
correlation function of the current operators S 3 defined 
in Eq. ifTTl . see also Refs. 

Let us now disregard interaction effects. In this case 
one should suppress fluctuations of ip in the end of the 
calculation by setting ip + {r) — eVr and ip~(r) — > 0. 
Then for the noise correlator one easily finds 



S 2 



(*e) 2 



Sip (ti)Sip (t 2 ) 



(12) 



and a similar expression is obtained for S 3 . Of interest 
is the irreducible part of the correlator £3 which reads 



S 3 = 



S 3 S 



(13) 



Sip (ti)Sip (t 2 )5ip (t 3 
= S 3 (t!,t 2 ,t 3 ) - (7(ii))«S(t 2 ,t 3 ) - (7(i 2 ))5(ti,t 3 ) 
-(/(t 3 ))5(t 1> t a ) + 2</(t 1 ))(/(t 2 )>(J(t 3 )). 

It follows immediately that in order to evaluate the third 
current cumulant in the absence of interactions it suffices 



to expand the exact effective action (J7J up to the third 
order in ip~ , 



iSo^] = iSV>[<p k ] + J 5< 2 )[^ ± ] +iS^[ip^ 



(14) 



keeping the full non-linearity in ip + in each of these terms. 
The first two terms of this expansion were evaluated in 
Ref. 0. Being combined with Eq. fL^I . the term S^ 2 ' [(^ ± ] 
allows to recover the well known expression for the shot 
noise spectrum 1 . Proceeding further with the expansion 
in ip~ for the term iS^ one finds 10 



dr(ip-(T)fip+( T ) 



(15) 



dn / dr 2 / dr 3 ip (n)ip {r 2 )<p (r 3 ) 



2nij 

3e 2 R J j j a 

Xf(T 2 ,Tl)f(T 3 ,T 2 )f(n,T 3 ), 

where 1/R = (2e 2 /h) J2 n T n is the scatterer conductance 
and 

f( , Tsin[(^+(r 2 )- y +(r 1 ))/2] 
/(T2 ' Tl) = sinh^T^-n)] ' (16) 

Let us substitute the above expressions into Eq. Ijl3|l 
and, after taking derivatives over ip~ , set ip + (r) = eVr 
and ip~ — > 0. This is sufficient provided the time dif- 
ferences \t\ — i 2 |, \t\ — t 3 \ exceed the charge relaxation 
time trc and provided eV <C ^)trc- The formula 1|13|) 
immediately yields 



S 3 = f3e 2 IS{h - t 2 )5(t 1 - t 3 ) 



(17) 



R 



f(t 2 -t 1 )f(t 3 -t2)f(h-t 3 ), 



where /(r) = Tsin(eV"r/2)/ sinh(7rTr). Performing the 
Fourier transformation 

5 3 K,W2)= / dr 1 dT2e i ^ +i ^S 3 (t 1> t 1 -n,t 1 -T2) 

we arrive at the final result 

S 3 = Pe 2 I -2 1 e 2 IF(v,w 1 ,w 2 ), (18) 

where 

sinh 3 (w/2) f 00 dtu 



F = 



(19) 



4« J-oo x(w)x(w-Wi)x(w + tc 2 ) 

Here we defined v — eV/T, w\, 2 — wi, 2 /2T and 

x(lu) = cosh 2 lu + sinh 2 (w/4). (20) 

Eqs. I|18|) - (|2U|I represent the main result of this paper. 
They fully describe the third cumulant of the current 
operator at voltages and frequencies smaller than both 

1/trc an d 1/td. 

Let us briefly analyze Eqs. Q18 )) -Q2U [I in various limits. 
For u>\ 2 = we recover the well known result 



F(v,0,0) = 1 + 3 



1 — (smhv/v) 
coshw — 1 



(21) 
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FIG. 1: The function F(v,w,Q) = F(v,w, -w), see Eq. fify . 



which in turn yields F — > 1 in the limit of large voltages 
v » 1. Eq. also holds for iu lj2 < t>. 
In the limit u<l one finds 



F(v < 1, io, 0) = F(u < 1, w, -w) = 
9 sinh w + sinh(3w) — 12w cosh w 2 



(22) 



48 sinh 5 ■ 



and, similarly, 



F(v<l,w,w)= (23) 
sinh(4w) + 4sinh(2w) — 8u; cosh(2w) — 4w 2 
128 sinh 5 tu cosh 3 iu 

These equations demonstrate that at large frequencies 
w 3> 1 the function F decays exponentially with w. From 
Eqs. (HD and (|53J) we find respectively F oc u 2 e" 2u 73 



Finally let us turn to the most interesting limit of low 
temperatures, in which case one always has v,Wx,w 2 3> 
1. Neglecting small corrections ~ 1/v, 1/w we obtain 



F(v,w 1 ,w 2 ) 



Wl2 



i{2\w 12 \<\v\ (24) 



F(v,w u w 2 ) = 0, i£2\w ia \>\v\. (25) 



Here the value W12 is defined differently depending on 
the sign of the product w\w 2 . For w\w 2 > we have 
w\ 2 = w\ + w 2 while in the opposite case wiw 2 < we 
define wi 2 — max \w 2 \]. We observe that in both 
cases the function F depends linearly on frequency and 
vanishes as soon as |wi2| exceeds \v\/2. 

At arbitrary values of v, w\ and w 2 the integral 119|) 
can be evaluated numerically. The corresponding result 
for the function F(v,w,0) is depicted in Fig. 1. The 
overall form of the function F(v,w,w) is similar but - 
as compared to F(v,w,0) - it demonstrates a somewhat 
faster decay with increasing frequency. 

In conclusion, we have presented a general approach 
which allows to describe statistics of current fluctuations 
in mesoscopic coherent conductors at arbitrary frequen- 
cies and in the presence of interactions. Restricting our- 
selves to the non-interacting case, we have analyzed fre- 
quency dispersion of the third cumulant of the current 
operator. This dispersion was found negligible only in 
the case of tunnel junctions, while in a general case it 
turns out to be important in the frequency range compa- 
rable to eV. Similar results are also expected for higher 
order cumulants of the current operator. 
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